Ordering in weakly coupled random singlet spin chains 
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The influence of bond randomness on long range magnetic ordering in the weakly coupled S = 1/2 
antiferromagnetic spin chain materials Cu(py)2(Cli_a,Br x )2 is studied by muon spin rotation and 
bulk measurements. Disorder is found to have a strong effect on the ordering temperature Tn, 
and an even stronger one on the saturation magnetization mo, but considerably more so in the 
effectively lower-dimensional Br-rich materials. The observed behavior is attributed to Random 
Singlet ground states of individual spin chains, but remains in contradiction with chain mean field 
theory [1] predictions. In this context, we discuss the possibility of a universal distribution of ordered 
moments in the weakly coupled Random Singlet chains model. 



The ground states of unfrustrated classical systems are 
typically robust with respect to weak Hamiltonian disor- 
der. A case in point is the classical Heisenberg antiferro- 
magnet (HAF). Randomizing the strength (but not the 
signs) of exchange interactions leaves the fully aligned 
Ncel ground state completely intact. In contrast, in 
quantum systems, arbitrarily weak disorder will modu- 
late the strengths of local quantum fluctuation and of- 
ten qualitatively reconstruct the ground state. The one- 
dimensional quantum S — 1/2 HAF is an extreme exam- 
ple. For uniform chains, the ground state is a Tomonaga- 
Luttinger spin liquid (TLSL) [2]. The introduction of 
arbitrary weak bond randomness gives rise to the so- 
called Random Singlet (RS) phase [3-6]. In the RS state, 
spin correlations are protected from localization effects 
by particle-hole symmetry [2, 7]. Nevertheless, the scal- 
ing laws [6, 8, 9], although universal and independent of 
the details of disorder, are markedly different from those 
of the TLSL. In experiments on real materials, one has 
to deal with quasi-one-dimensional (quasi-lD) spin sys- 
tems. A divergent correlation length in individual chains 
ensures 3D long range order at Tn > for arbitrary weak 
inter-chain interactions J', in both the disorder- free [10] 
and disordered cases [1, 6]. An intriguing question is to 
what extent the peculiarities of the RS state in d = 1 
translate into unusual features of the ordered phase in 
d = 3. 

To date, RS-forming bond randomness in quasi-lD 
magnets has received considerably less attention than 
dilution- type disorder due to spin substitution [11-13]. 
For the former, the existing predictions are derived from 
chain-Mean Field (chain-MF) theory [1]. The main re- 
sult is that random bonds tend to increase both Tn and 
the ordered moment m at T — > 0. This effect is related 
to an abundance of very loosely coupled and almost free 
spins in the RS state of isolated chains [3, 5]. In coupled 
chains, for weak J', one gets Tn oc J' mo- In contrast, for 
disorder-free case, the in-chain interactions enter the rela- 
tion explicitly: Tn oc JrriQ [14]. On the experimental side, 
the challenge is to measure the very small sublattice mag- 



netization that arise in the weak-coupling regime, where 
this theory may be expected to apply. In the present 
work we overcome this difficulty by employing the sen- 
sitive muon spin rotation (/iSR) technique, which has 
emerged as a tool of choice for the study of quantum mag- 
netism [15]. We directly measure the relative variations 
of mo and Tn in the prototypical bond-disordered quasi- 
lD S = 1/2 HAF systems Cu(py) 2 (Cli_ x Br x ) 2 . We find 
that both quantities are reduced by disorder and that mo 
is far more sensitive than the transition temperature, in 
apparent contradiction with chain-MF predictions. 

Our target compounds are derivatives of Cu(py)2Cl2, 
one of the first known and extensively studied 5 = 1/2 
spin chain materials [16]. Single crystals samples with 
varying Br content x are straightforward to grow from 
solution by slow evaporation. In Cu(py)2Ci2 (space 
group P2 1 /n, a = 16.967A, b = 8.5596A, c = 3.8479A, 
j3 = 91.98°) the chains are formed by magnetic S = 1/2 
Cu 2+ ions linked by superexchange bonds via the halo- 
gen sites. The temperature dependence of magnetic sus- 
ceptibility (Fig. la) shows a broad Bonner-Fischcr (BF) 
maximum [17, 18], characteristic of a quantum 5=1/2 
chain with an AF exchange constant J x —o — 2.35 meV. 
The chain-structure of Cu(py)2Br 2 is quite similar ( 
P2 1 /n, a = 8.424A, b = 17.599A, c = 4.0504A, = 
97.12° [19]), as are the measured magnetic susceptibil- 
ity curves. However, the in-chain exchange constant is 
larger: J x= \ = 4.58 meV. Weak inter-chain interactions 
lead to 3D ordering in both materials, at Tn = 1.15-ftf K 
[20] and T N = 0.72 K for Cu(py) 2 Cl 2 and Cu(py) 2 Br 2 , 
respectively. The transitions are marked by well-defined 
lambda-anomalies in the measured temperature depen- 
dence of specific heat C(T), as shown in Fig. lb [33] 
Knowing T N and J, allows us to estimate the effective 
inter-chain coupling constants [21]: J' x=a = 0.05 meV 
and J' x= i = 0.03 meV for the two materials, corre- 
spondingly. The bromide is clearly a much more one- 
dimensional system. This is also evident from its much 
weaker lambda-anomaly. It is useful to estimate the or- 
dered moment at T — ► 0. Based on chain-MF results, for 
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FIG. 1: (Color online) Bulk properties of Cu(py)2(Cli_ a: Br I )2 
samples, (a) Measured temperature dependence of magnetic 
susceptibility for a field applied along the chain axis (sym- 
bols) , and fit of the theoretical curve for the uniform quantum 
S = 1/2 HAF chain [17, 18] (solid line). From the buttom 
up, the data are offset along the y axis by 0, 1, 2, 3, 4.2, and 
5.8 10 -3 emu/mol, respectively, (b) Measured temperature 
dependence of specific heat (symbols). The vertical offsets 
are 0, 0.16, 0.27, 0.55, 0.75 and 1 J/mol K, respectively. 



Cu(py) 2 Cl 2 and Cu(py) 2 Br 2 we get m , x=0 = 0.15 ^ B 
and mo tX= i = 0.08 /iB, respectively. 

As determined by single crystal X-ray diffraction and 
chemical elemental analysis, the structures of Cu(py) 2 Cl 2 
and Cu(py) 2 Br 2 are stable with respect to chemical sub- 
stitution of CI for Br and vice versa, for x < 0.4 and 
x > 0.6, respectively. Due to a difference in CI and Br 
ionic radii, such disorder randomizes the Cu-halogen-Cu 
bond angles and thereby the magnitude of in-chain ex- 
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FIG. 2: (Color online) Composition dependence of the in- 
chain exchange constant J (from fits to x(T) ), the ordering 
temperature Tn (from calorimetry) and the low-temperature 
ordered moment mo (chain-MF estimates and /^-SR measure- 
ments). 



change interactions, whilst leaving the spins themselves 
unaffected. This strategy of creating bond-disordered 
systems has been previous successfully applied in other 
materials such as IPA-Cu(Cli_ 2; Br a .)3 [22], piperazinium- 
Cu^Cl^Br^e [23], H s C 4 S0 2 -Cu 2 (Cl 1 _ ;c Br ;c ) 4 [24] and 
NiCl 2 -4SC(NH 2 ) 2 [25]. Due to the slightly differ- 
ent structures of the parent compounds, we are actu- 
ally dealing with two tunable random-bond spin chain 
materials, on the Cl-rich and Br-rich ends of the 
Cu(py) 2 (Cli_ a .Br :E ) 2 line, respectively. We will reserve 
the formula Cu(py) 2 (Cli_ x Br2,) 2 for Cl-rich compounds 
(x < 0.5) and use Cu(py) 2 (Bri_ y Cl a ) 2 to denote materi- 
als on the Br end (y < 0.5). 

The bulk properties of the halogen-disordered sam- 
ples resemble those of the corresponding disorder-free 
systems. Typical measured magnetic susceptibility data 
are plotted symbols in Fig. 1. The derived average in- 
chain exchange constant J steadily increases with in- 
creasing Br content. A 3D magnetic ordering transi- 
tion is observed at low temperatures in all composition 
studied(Fig. 1, lower panel). Tn decreases with Br con- 
tent on the Cl-rich end and with CI concentration in Br- 
rich samples (Fig. 2). The variation is more pronounced 



in Cu(py) 2 (Br!_ a Cl y ) 2 . Moreover, the corresponding 
lambda-anomaly weakens and slightly broadens with in- 
creasing CI content in Cu(py) 2 (Bri_j,Cl a ) 2 , while it re- 
mains almost unchanged in Cu(py) 2 (Cli_ x Br a; )2. The 
measured transition temperature and average exchange 
constant are plotted versus composition in Fig. 2. 

The main focus of the present study is on ^i-SR mea- 
surements. This technique has been instrumental in the 
study of very small static[26, 27] and dynamic[28, 29] mo- 
ments in spin chain systems. It probes the local magnetic 
fields at the stopping sites of muons implanted into the 
sample [30]. These, in turn, are expected to be propor- 
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FIG. 3: (Color online) Muon spin asymmetry measured vs. 
time in several Cu(py)2(Cli- 2 ,Br :E )2 and Cu(py)2(Bri_j / Clj / )2 
samples at T ~ 20 mK (solid symbols) and just above Tn 
(open symbols). The plots for y = 0.05 and 0.1 and x = 
0.05, are offset along the y axis by 0, 0.1, 0.2, 0.3 and 0.43, 
respectively. The solid lines are fits to the data as described 
in the text. 

tional to the static ordered moment. We performed mea- 
surements on powder samples of Cu(py)2(Cli_ a; Br a ;)2 
and Cu(py) 2 (Br 1 _j / Cly)2 with x — 0,0.05 and y = 0, 
0.05 and 0.1 at the LTF spectrometer at the S/zS muon 
source at Paul Scherrer Institut. Typical muon spin re- 
laxation/rotation curves measured in zero applied field 
(ZF) are shown in Fig. 3. Several distinct time scales 
are apparent. In all samples, a very rapid decay at 
short times can be attributed to muonium formation with 
the organic ligand [31]. This contribution appears to be 
temperature-independent and is not directly relevant to 
the physics discussed here. At the largest time scales, a 
slowly decaying tail is due to the non-precessing muon 
spin components parallel to the local field, and to muons 
stopping outside the sample [30]. 

As expected, in the paramagnetic phases of all samples 
one only observes a slow decay of muon polarization, with 
no oscillatory behavior. Representative data collected at 
T > Tn are shown in open symbols in Fig. 3. They can 
be modeled with exponential decay processes, as shown 
in dashed lines. We conclude that above the ordering 
transition, all materials studied behave very similarly, 
despite the different levels of disorder. 

Clear differences emerge at low temperatures, in the 
magnetically ordered phase (Fig. 3, solid symbols). Al- 
though the relevant time scales turned out to be strongly 
dependent on composition, it is possible to provide a 




FIG. 4: (Color online) Magnetic order parameter (charac- 
teristic field at the principal muon stopping site) plotted 
against temperature in units of the average in-chain exchange 
constant J in Cu(py)2 (Qi-ajBr^ and Cu(py)2(Bri_ !/ Clj / )2 
(symbols). Within each family of materials, the relative mag- 
nitudes of the order parameter can be directly compared. 
Lines are guides to the eye. 

common description of the low-temperature /i-SR spec- 
tra in all samples. In addition to the fast and and slow 
background contributions described above, for T < Tn, 
the main effect is the appearance of spontaneous muon 
spin precession. An application of a small longitudinal 
field parallel to the muon spin direction recovered all po- 
larization apart from the muonium contribution. This 
observation shows that the observed spin relaxation and 
rotation are due to static internal fields, which we at- 
tribute to ordered Cu 2+ moments. We chose to model 
these processes as a sum of a damped Bessel function and 
an exponential term. The spectra are then described as: 

A(t) = A 1 J a (ujt + <f>)exp(-\ 1 t) + A 2 cxp(-\2t) + 

+ Afest exp(-Afastt) + Ataiiexp(-A ta il*)- (1) 

The first term on the RHS represents the precession 
of muons stopped at the most probable sites [28]. The 
choice of a Bessel function, typically used to describe 
incommensurate structures [30] is in our case purely em- 
pirical. Nevertheless, its use is justified by preliminary 
neutron diffraction evidence [34] that the magnetic struc- 
ture is actually helimagnetic in the b direction, perpen- 
dicular to the chains. Its exponential envelope reflects 
a narrow Lorentzian distribution of local fields and also 
takes into account the depolarization by nuclear spins. 
The characteristic muon spin precession frequency uj is 
directly proportional to the magnitude of the static mag- 
netic field uj = 7-B, with 7 = 85.16 krad s _1 G _1 , is thus 
our primary measure of the static magnetic order. 

The second term in Eq. 1 empirically describes the 
multitude of other stopping sites sensitive to the static 
magnetic order. For each sample, the site occupancies 



4 



A\ and A2, as well as parameters for "fast" and "tail" 
contributions, were determined in global fits to the data 
collected in the entire temperature range. The parame- 
ter 4> determines the functional shape of the oscillatory 
term. At each temperature, it was globally applied to all 
samples, separately on the Br-rich and Cl-rich ends, to 
allow a direct comparison of the precession frequencies u) 
within each family of materials. For each sample, the pa- 
rameters oj, Ai and A 2 were refined at each temperature 
separately The model provides excellent fits to all data 
collected for the ordered state in all samples. Typical fits 
are shown in heavy solid lines in Fig. 3. The values of 
all fit parameters for all samples and their temperature 
dependencies are deposited as supplementary material. 

The main result of our analysis is the temperature de- 
pendence of the order parameter B, plotted in Fig. 4. 
For each composition, the temperature axis has been 
normalized to the value of T c derived from calorimetric 
measurements. Three conclusions can be immediately 
drawn. First, in apparent contradiction with chain-MF 
theory, both T N and the saturation magnetization to 
are reduced in disordered samples, on both composition 
ends. Second, the Br-rich materials, which are more one- 
dimensional in nature, are considerably more affected by 
disorder that the more 3-dimensional Cl-rich systems. In 
Cu(py) 2 (Br 1 _ J) Clj,)2, the oscillatory term is dramatically 
slowed already at 5% substitution, and is almost com- 
pletely overrun by relaxation effects at 10% CI. Third, 
for Br-rich samples the effect of bond disorder on too is 
much more drastic than on Tn. This is made particu- 
larly clear by the too vs. composition plot in Fig. 2. It 
is based on our initial estimates for too in the disorder- 
free materials, and on the assumption that within each 
of the two series of materials too is proportional to the 
Larmor field B extrapolated to zero temperature. Note 
that for T — > the ratio of observed fields for the two 
disorder-free samples is actually very close to the ratio of 
estimated ordered moments, so the data in Fig. 4 can in 
practice be compared across all compositions. 

Before we discuss these results, a few words need to be 
said on the the nature of the bond disorder. That disor- 
der effects are not apparent in bulk properties above Tn is 
not at all surprising. Recent numerical simulations and 
experimental studies of the RS material BaCu2SiGe07 
[32] have illustrated that the effect of even very large 
randomness of exchange constants may be very modest, 
as far as bulk properties are concerned. In fact, de- 
viant behavior of susceptibility and magnetization curves 
may emerge only at temperature and field scales that are 
over an order of magnitude smaller than the width of 
the random bond distribution. Even if we assume that 
in Cu(py)2(Bri_j / Clj,)2 we arc dealing with an extreme 
and clearly exaggerated bimodal distribution of exchange 
constants, a 10% substitution will produce a standard de- 
viation of only about 0.7 meV. In this conservative esti- 
mate, RS behavior may be expected to affect the bulk 



properties only at temperatures well below ~ 0.5 K, 
where we will have already entered the ordered phase. 
Here it should also be noted that the effect of chemical 
disorder on mier-chain interactions will have a minimal 
effect on 3D ordering. It will be averaged out by the 
very large correlation length in individual chains at low 
temperatures. 

The effects of disorder and the drastic suppression 
of to relative to T N are most prominent in the Br- 
rich materials, which are more one-dimensional. Here 
the RS-related physics is more likely to play a role. 
RS properties emerge only below a certain energy scale 
that is non-universal and depends on the initial dis- 
tribution of exchange constants [3-6]. In contrast to 
Cu(py) 2 (Br!_ y Cl y ) 2 , in Cu(py) 2 (Cli_ a: Br a; )2 this energy 
scale may be much lower than J', making the disorder ir- 
relevant in the ordered phase. The discrepancy between 
our findings and the chain-MF theory of Ref. [1] remains 
to be explained. One tantalizing possibility is that the 
chain-MF approach may be in principle be inapplicable 
to the RS phase [35] . The latter features an abundance of 
weakly-dimerized spin degrees of freedom for which quan- 
tum correlations with similar objects in adjacent chains 
simply can not be ignored. 

A final note concerns the homogeneity of the static or- 
dered moment. The low energy physics of the RS phase 
is exactly that of non-interacting random dimers with 
a universal probability distribution of dimcr strength. 
These singlets will be partially polarized by the mean 
exchange field. The degree of polarization will be deter- 
mined by the strength of the dimers. As a result, we 
expect a universal probability and spatial distribution of 
static ordered moments P{m) in the T — > limit. Even 
though muon spectroscopy could, in principle, measure 
this distribution directly, in Cu(py) 2 (Cli_ a; Br 2 ;)2 and 
Cu(py) 2 (Br 1 _ y Clj / )2 such an experiment appears very 
challenging. We believe that this elegant idea of a univer- 
sal distribution of ordered moments in weakly coupled RS 
chains deserves more theoretical and experimental atten- 
tion, using new materials for ^-SR and using alternative 
techniques such as NMR. 
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